The general relativistic geometry of the Navarro-Prenk- White model 
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We derive the space time geometry associated with the Navarro Frenk White dark matter galactic 
halo model. We discuss several properties of such a spacetime, with particular attention to the 
corresponding Newtonian limit, stressing the qualitative and quantitative nature of the differences 
between the relativistic and the Newtonian description. We also discuss on the characteristics of the 
possible stress energy tensors which could produce such a geometry, via the Einstein's equations. 
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The geometry generated by the galactic halo, is generally thought to be "almost flat". This assumption is based, 
first, on the fact that the galactic dark halo has very low density, at most a few order of magnitude above the critical 
one. Second, the velocities involved are small compared with the speed of light, and third, the dust treatment gives a 
description of the dynamics which is in good agrement with the observation. Thus, validates the Newtonian physics 
\& ' as an adequate to treat the dark halo. 

These very same arguments are used for studying the Solar system, where also the geometry is taken as " almost flat" . 
Nevertheless, the general relativistic treatment of the Solar system has made possible to give important corrections 
to the Newtonian one and, furthermore, in the general relativistic treatment is precisely this "almost not flatness", 
which explains the motion of the planets! 

We consider that to count with a general relativistic version of the galactic dark matter halo allows one to make 
, more accurate analysis on the dynamics of the objects, including the study on gravitational lenses, to mention an 
C*~) ' application. 

In the present work, we describe how the observations can be related with part of the geometry, then propose and 
expression for the complete geometry associated with the Navarro-Frenk- White, NFW, model. Next describe the 
properties of such a geometry and explain why the Newtonian description works so well. Armed with the geometry, 
we discuss on the type of matter-energy which generates the geometry, that is, the nature of dark matter, a point 
Oh' where the Newtonian analysis remains mute. 

In a previous series of works, Matos et al.0, Guzman et. al. Q, we discussed the possibility of determining the 
geometry of the space-time, and then constraining the type of matter-energy which generate such geometry, based on 
observational data. In particular, we addressed the problem of making those determinations based on the observed 
C$ • profile of the tangential velocities of objects orbiting galaxies. 

In a nutshell: Given the fact that the dark halo in the galaxies seems to be spherical and at rest, at least in the 
average, we consider a general spherically symmetric static space-time, see Eq.Q below, and were able to determine, 
on purely geometrical ground, an expression for the tangential velocity of objects moving in circular stable geodesies 
in terms of the metric coefficients, which turned out to be a very simple one, Eq.@. We then took a sort of inverse 
point of view. Instead of consider such an equation as an expression for the velocity, we took it as an expression for 
the metric coefficient, given the fact that what is being observed is the velocity profile, thus being able to determine 
part of the geometry based only on observational data, Eq. (f 1 01) . 

In the present work we present such a program applied to the NFW model, Navarro et al. Q , which has proved to 
have a remarkable predicted power and agrees very well with observations, particularly with those outside the central 
galactic region. 

In what follows we reproduce the main steps on the reasoning leading to the conditions which the tangential velocity 
of circular orbits impose on the metric coefficients for the spherically symmetric static case. 
We begin with the general line element for such geometry: 

ds 2 = -a 2 (r)c 2 dt 2 + - / — v +r 2 dQ 2 , (1) 
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where c is the speed of light and G the gravitational constant, and dVl 2 = dd 2 + sin 9 dip 2 is the solid angle element. 
From the corresponding Lagrangian for a test particle in this space, 2C = (jr) > where r stands for the proper time, 
we obtain that the energy, E — a 2 (r)c 2 i, the (^-momentum L v = r 2 sin 2 #(p, and the total angular momentum, 

L 2 = Lg 2 + (-^g) 2 , with Lg = r 2 9, where dot stands for derivative with respect to the proper time, are conserved 
quantities along the motion. Notice that we can write the total angular momentum in terms of the solid angle as: 
L 2 =r 2 tl 2 . 

With this information, the fact that the four- velocity, = is normalized, u v u v = — 1, translates into a radial 
motion equation: 



with the potential V(r) given by 



+ V(r) = 0, (2) 



Restricting the radial motion to circular stable orbits, implies imposing the conditions, r = 0, for circular orbits, 
and — 0, so that it describes an extremum of the motion, and > 0, so that the extremum is a minimum. These 
three conditions guarantee that the motion will be circular and stable. They also imply the following expressions for 
the energy and total momentum of the particles in such orbits: 

E* = ,f« 3 < r > , (4) 
a[r) — r a(r) r 

a[r) — r a{r) r 

where a subindex r stands for derivative with respect to r. 

On the other hand, we can rewrite the line element for this geometry, Eq.Q, in terms of the modulus of the 
spatial velocity, normalized with the speed of light, measured by an inertial observer far from the source, as ds 2 — 
—dt 2 (l — v 2 ), where 

v = c 2 a 2( r ) I " 9 a m(t\ ~^ r ( ~ ) ) ■ (^) 

This last equation implies that the modulus of the angular velocity, which is the tangential velocity for the case of 
circular orbits, is defined as: 
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c 2 a 2 (r) \ dt J c 2 a 2 (r) \dt J 
thus, in terms of the conserved quantities, the angular velocity takes the form: 



n 2 , (7) 



2 _ c 2 a 2 (r) L 2 
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Using the expression derived for these conserved quantities, Eq.|J5J, we obtain that the tangential velocity can be 
expressed in terms of the metric coefficient a as: 



V = (9) 
a(r) 

This last equation allows us to determine the metric coefficient a(r) in terms of the observed velocity profile: 

f v 2 (r) 

a(r) = cxp / JSli. dr. (10) 

This is the key equation of the reasoning: To use the observations, w tg (r), in order to partially determine the 
geometry of the surrounding spacetime. 
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Now we join these results with the Navarro- Frenk- White model. This model predicts the density profile [!|: 



Po 

^(i + T:) 2 



PNFW = r T-X I r_\2 > (^) 



where po = p cr n 5 C , r s is a scale radius, S c is a characteristic (dimensionless) density, and p cr i t = | ^ G is the critical 
density for closure. The mass function, Mnfw(^) = 4-7T J r 2 pNFW rf'') with the integration constant chosen so that 
Mnfw(?" = 0) = 0, takes the form: 



M NFW (r) = Anri Po \ hi [! + -)- Y^VJ ■ ( 12 ) 
This implies, equating the gravitational force with the centrifugal one, the following profile for the tangential velocity: 



'Mi + £) i 



tSNFW 



(13) 



where «q = 4 tt G po rl is three times the tangential velocity of particles at the limb of a sphere of radius r s and a 
mass defined by the critical density times the S c factor inside that volume. 

These expressions are directly predicted by the NFW model and, with exception of the central parts of galaxies, it 
has been successfully compared with the actual observations 0. 

Using the expression derived for the tangential velocity within the NFW model, Eq. (|13|) . in the expression we 
obtained above, Eq. @, we obtain a remarkable simple expression for the gu coefficient: 



a 2 {r) 



(14) 



where we have set the integration constant such that a goes to one for large radii, and we have normalized with the 
speed of light the NFW velocity. There are several noticeable features in this last expression. First, it is regular 
everywhere. The divergency problem that the NFW-density has in the central region, is not reflected in the metric 
coefficient: 



lim a 2 = e 



, lim a 

r — >oo 



1. 



(15) 



Actually, the a-function goes to one for large radii, as can be seen in the figure, ^ recovering and validating the 
Newtonian assumption in that region. 
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FIG. 1: The g t t metric coefficient, for values ^IJ- = 10 4 , r s = lOfcpc. 
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We can go on and work with the other metric coefficient, g rr = (I — 2G ^ r J- r ^ First of all, recall that the gu 
coefficient was determined based on the analysis presented above. For the g rr coefficient, the unknown function M(r) 
can not be directly identified with the mass function. Strictly speaking, the mass is defined for asymptotically flat 
spacetimes with several characteristics, but we do not go into that discussion in this work. We can say that this is 
an approximation. On the other hand, however, as will be shown below, the \ component of the Einstein equations, 
implies a consistency between the fact of taking this function as the mass and the defined density, so that we have 
some grounds to consider it like the mass. Thus, let us be bold and take the function M(r) as the mass function of 
the NFW model. The line element takes the form: 

ds 2 = - (l + -) 2 W c 2 dt 2 + , ^ + r * ^ ( 16 ) 

H 1+ t-) i_ x 

77 1 + 77 

In this way, we obtain the geometry associated with the Navarro-Frenk- White model. 

The line element has two free parameters, namely Do, the characteristic speed, and r s , the characteristic radius. 
With respect to the g rr metric coefficient, it is also regular everywhere, for positive radius, we see its behavior in 
figure Fig.©, and it has the following limits: 

lim g rr = lim g rr = 1. (17) 

r — >0 r — >oo 
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FIG. 2: The g t t metric coefficient, for values = 10 4 , r s = Wkpc. 



Thus, the line element given by Ea. <|16f) is regular at all points and far from the source takes the form of the flat 
spacetime. 

However, even though the line element does not show any divergence, we expect to have one, directly inherited 
from the divergence of the NFW density at the center. Actually, when one derives the scalar of curvature, R, it is 
seen that: 

Rr^o ~ 2^L - f 2 - 4) + OW, (18) 

showing a clear divergence at the center, as expected. This implies that we do not have a solution describing a 
spacetime, as long as there is a naked singularity, the divergence is not covered by any horizon. What we have is a 
geometry which implies a dynamic similar to the one described by the NFW model, so we claim that this geometry 
will be related with the exterior part of a complete spacetime with the observed dynamics. 

The metric coefficients differ from the unity by very tiny amounts, validating the flat space approximation, that is, 
the Newtonian analysis. However, within the General relativistic theory, these tiny amounts of non-flat geometry, are 
the ones responsible for the observed dynamics! 
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Finally, we can construct the Einstein tensor and obtain some conclusions about the type of matter-energy which 
produce the space time given by Ea. (|l(jf) . The Einstein tensor gives three non-zero independent components: 

G\ = _2-^-L (19) 

c 2 r 2 s xu 2 

v 4 1 

G r = ~ A ^T2 -ITS [ ulnu (u^u-2x) + x 2 ] , (20) 



G 6 e = (ulnu- 



2u\nu [%ulnu-2x (u+%\) + x 2 (^7a; + 4 + 2^ 



(21) 



where we have defined x = — , and u = 1 + — . 

Being aware of the caveats on promoting this geometry to a spacetime, still we can say something about the type 
of matter which could produce such a geometry, by means of the Einstein equations: 

G» = ^-T», (22) 
c 4 

where G stands for the gravitational constant and T£ describes the tensor of distribution of the matter-energy in the 
spacetime. 

It is common to identify the T t * component of the matter-energy tensor with the density of matter, and energy, 
present in the spacetime, p, that is T t * = — c 2 p. Thus, using the corresponding expression for the Einstein tensor, 
Eci. (|19fl . we obtain the same expression relating the mass and the density as that obtained in the Newtonian theory, 
Ea. (|12fl . This fact gives support to the interpretation of the function M(r) in the line element, Ea. l|16|) . as the mass 
of the system. These are good news and give consistency to the treatment presented here. 

About the other components of the matter-energy tensor, by means of the Einstein equations, we can conclude that 
such matter-energy tensor can not be dust or even perfect fluid. The reason for these conclusions are clear. Again, if 
we take the matter-energy tensor to be a perfect fluid one 

T$= (pc 2 +p)u^u l/ +p8^ (23) 

with u M the four velocity. For the spherically symmetric case we are dealing with, the matter-energy tensor takes 
the form Tjf — diag (— pc 2 ,p,p,p). But, from the Einstein tensor, we see that clearly G r rl and G 9 g are non zero and 
different. Thus, the matter-energy curving the spacetime to form the NFW geometry can not be dust or perfect fluid. 

If we insist on a dark fluid type for the matter-energy, it has to be an anisotropic one, with two different pressures, 
a radial one, p r , and a tangential one, p± which, by the Einstein equations imply 

9 9 

p r = -p 2 2nG-^ r - [ulnu (u lnu - 2 x) + x 2 ] (24) 



p± = —p 2 8irG-^r- (ulnu-i) 
x- 1 
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c 2 \ c 2 J J V c 2 



(25) 



It is important to notice that, as the density, p is very small, these pressures are even smaller, and tend very quickly 
to zero. Thus, again validate the Newtonian treatment taken the fluid as dust, in an exterior region. Nevertheless, 
recall that these analysis is on the track of determining the actual nature of the dark matter, pointing to the physical 
properties it must have. 

In this way, the geometry associated with the NFW model that we have derived, allows us to explain the validity 
of the Newtonian treatment, pointing to its limits, to discuss on the nature of the dark matter, and opens the road 
to perform several dynamical studies within the general relativistic theory. 
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